Anderson localization due to random magnetic field in two dimensions 



OO 
On 
On 

> 
o 

(N 



Akira Furusaki 

Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 
and Department of Physics, Stanford University, Stanford, California 94305 

(February 7, 2008) 

Results of large-scale numerical simulations are reported on the Anderson localization in a two- 
dimensional square lattice tight-binding model with random flux. Localization lengths, fluctuations 
of the conductance, and the density of states are computed for quasi-one-dimensional geometry. 
Numerical results indicate that the model exhibits the same critical behavior as the one studied 
by Gade and Wegner. It is argued that all the states except a zero-energy state are localized and 
the density of states has a singularity in the center of the band. The energy scale below which the 
density of states increases is found to be extremely small (< 10 -2 ). 
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It is general wisdom that noninteracting electrons are 
localized in two-dimensional (2D) disordered systems [Q. 
There are, however, some well-known exceptions to this 
rule. These include electrons having strong spin-orbit 
coupling H and integer quantum Hall systems |!|. Re- 
cent studies have shown that 2D Dirac fermions with ran- 
dom gauge field offer another exception to the rule . 
For a model of 2D nonrelativistic fermions subjected to 
random magnetic field with zero mean, the existence of 
delocalized states has been a subject of debate. 

The random flux model, in which static magnetic field 
is randomly distributed with zero mean, got much atten- 
tion recently in connection with the gauge field theory of 
high-r c superconductivity |J and the composite-fermion 
theory of the half-filled Landau level Q. It has been 
controversial, however, whether this model has a delo- 
calized state ||. On the one hand, several numerical 
and analytical studies concluded that all the states are 
localized and belong to the unitary class of the scaling 
theory p|-p^]. On the other hand, a different conclusion 
that there are delocalized states near the center of the 
band was reached by other people ]f5|-p2[. One source 
of the discrepancy in numerical works is the extremely 
large localization length near the band center, making 
it difficult to decide whether or not states are localized 
from numerical data of finite-size systems. 

In this paper I present various numerical results ob- 
tained through the largest numerical simulations per- 
formed so far for the square lattice tight-binding model 
subjected to random flux with zero mean. The results 
indicate that a state at the band center (E = 0) is not 
localized. This is reminiscent of the integer quantum Hall 
system. There is, however, a crucial difference: the den- 
sity of states (DOS) is found to be divergent at E = 
in the random flux case. This behavior is similar to the 
ID and 2D random hopping models 



23 24 1 , and a crucial 
role is played by a special particle-hole symmetry relat- 
ing a state of energy E with a state of energy — E. The 
random flux model is argued to be in the same univer- 



sality class as a model studied by Gade and Wegner [B5| . 
Although this was already anticipated in Ji],^o| , this Let- 
ter reports for the first time that the random flux model 
shares a hallmark of the Gade- Wegner model, i.e., the 
divergence of the DOS at E = 0. 

The Hamiltonian of the tight-binding model is 



C j+l,k C 3,k 



" C j,k+l C 3, k 
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where Cj t k is annihilation operator of a fermion on site 
(J, k). The random magnetic flux is introduced through 
the random Peierls phase 6j t k in the hopping matrix el- 
ement. The magnetic flux <j)j ^ = Oj^ — Qj-i,k takes a 
random number in — up < (f>j^ < irp with a uniform dis- 
tribution. The parameter p is set to be 1, except in Fig. ^|. 
Numerical calculations are done for samples that have 
quasi-lD geometry of width M in the y direction and of 
length L in the x direction [M <C L). Periodic boundary 
condition is imposed in the y direction (c^m+i = c j,i): 
whereas open boundary conditions are assumed in the x 
direction for most of the calculation. For even M the 
lattice can be divided into A and B sublattices. For each 
eigenfunction ipE with energy E, changing sign of ipE on 
every site of, say, the A sublattice yields a new eigen- 
function tp—E with energy —E [^o| , ^6| . This symmetry 
relating the ±E states holds for each disorder configura- 
tion. For odd M, however, the particle-hole symmetry is 
absent under the periodic boundary condition. 

The localization length is calculated from the ex- 
ponential decay of the retarded Green's function ob- 
tained by using the standard recursive algorithm p7| : 
(In \\G r E (l,k;L, k') ||) ~ -L/X M , where ||G|| and ( ) de- 
note norm of G and ensemble average, respectively. Fig- 
ure Q shows the quasi- ID localization length Am normal- 
ized by M as functions of M and E. The typical length 
of quasi-lD samples used in the calculation is 3 x 10 5 , 
4 x 10 5 , and 8 x 10 5 for M = 32, 64, and 128, respec- 
tively. Furthermore, ensemble average is taken, typically, 
over 70 (20) samples for M < 64 (Af = 128) to reduce 
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the statistical error. The quality of the numerical data 
is therefore greatly improved from the earlier numerical 
results [p|, ^3| , ^6|JT^j20| ] . Clearly the states near the band 
edges (\E\ > 3.0) are localized [Fig. ||(a)]. Figure |](b) 
shows Xm/M decreases as M increases, suggesting that 
the states with \E\ > 0.1 are all localized in the 2D limit. 
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FIG. 2. Scaling curve obtained from the data for 
-3.5 < E < -2.55. For -3.0 < E < -2.55 only the data 
of M > 32 are used. Inset: Localization length versus energy. 




FIG. 1. Localization length for quasi-lD geometry calcu- 
lated for M = 4, 8, 16, 32, 64, 128, and 256. Additional 
data of M = 15 and 31 are shown for E = —0.4 and —1.0. 
In Fig. |l|(b) the statistical error for the data at M = 128 is 
about the same size as the symbols, whereas for smaller M 
the error bar is much smaller than the symbols. The kink at 
small M in the E = —0.1 data should be finite-size effects. 

The localization lengths of the quasi-lD wires are ex- 
pected to satisfy the one-parameter scaling Xm/M — 
f(£/M), where £ is the localization length in 2D. The 
scaling indeed holds as shown in Fig. || pq] . The scal- 
ing curve quantitatively agrees with the earlier results of 
Refs. H and JhJ. The agreement with the latter work 
is somewhat surprising in that a network model is used 
in which is an effective model in the semiclassical 
limit. The 2D localization length £ grows exponentially 
and reaches 10 6 lattice spacings at E = —2.55; see inset. 



FIG. 3. Localization length at E = as functions of M 
and the on-site disorder w. The filled symbols represent data 
for M = 7, 15, 31, 63, and 127. The open symbols are data 
for M = 8, 16, 32, 64, and 128. 

Figure | shows Xm/M versus M at E = 
There is a striking even-odd effect in the E = data, 
as noticed earlier in Refs. [^0|,^0). A new finding here 
is that Xm /M\ e _ q stays almost constant for odd M 
while it gradually increases for even M, suggesting that 
Xm/M — > const(> 0) as M — > oo. This would mean 
that ^.e=o is a critical or multifractal wave function, 
as suggested by Miller and Wang |2(J. By contrast, at 
\E\ =0.1, there is little even-odd oscillations [Fig. 0(b)], 
and Xm/M is a decreasing function of M. The impor- 
tance of the particle-hole symmetry can be also seen by 
examining the effects of on-site disorder, which breaks the 
symmetry. The on-site disorder is introduced by adding 
a term k effect k Cj.k to H, where Cj t k are taken to be 
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randomly distributed in the interval [— w/2, w/2]. Fig- 
ure H clearly shows that in the presence of the on-site 
disorder Am's of even and odd M's merge together and 
decrease for M > M c . The crossover width is M c w 64 
(300) for w — 0.4 (0.05) and would diverge as w — > 0. 
These results strongly suggest that the wave function is 
critical in the 2D limit only when E = w = 0, and that 
small | E | or w is sufficient to change it to a localized one. 

The states away from the band center belong to the 
unitary class. This can be verified by calculating fluc- 
tuations of two- terminal conductance as a function of L. 
For this purpose, perfect leads are attached to both ends 
of quasi-lD wires, and the transmission matrix t is calcu- 
lated from the Green's function G r E . The dimensionless 
conductance g is then obtained from the Landauer for- 
mula, g = Tr(tt^). Figure [i] shows vaxg — (g 2 ) — (g) 2 for 
M = 32, averaged over 2 x 10 4 samples. For \E\ = 0.1, 
varg is calculated for p — 1 and 0.2 without the on-site 
disorder. Almost identical varg versus L/Xm curves are 
obtained for \E\ = 1.0 and 0.02 as well. A thin line in 
Fig. ^ shows var g of the unitary ensemble calculated in 
the limit M <C L by Mirlin et al. [flj using the super- 
symmetric a model approach. Notice that, except for the 
peaks at L < 0.5Am, the numerical results of \E\ = 0.1 
are indistinguishable from the thin line (unitary ensem- 
ble). The discrepancy occurs only for L < M, where 
the samples are no longer quasi one-dimensional. The 
numerical curve of p = 0.2 is closer to the analytic re- 
sult because M/Xm\ p =o.2 <SC M / Xm\ p =i- These results 
clearly show that for \E\ > 0.1 and p > 0.2 the wave 
functions belong to the unitary class. 

The variance of g has a different L/Ajvi-dependence at 
E = for even M; see inset of Fig. f|. Without the 
onsite-disorder, for each L/Xm, varg of E — is larger 
than var g of E ^ j32|. This clearly shows that, when 
W = 0, the zero-energy state does not belong to the uni- 
tary class. The on-site disorder, however, drives ipE=o 
back to the unitary class, as shown by the long-dashed 
line (p = 1 & w = 0.2) in the inset. These observations 
are consistent with Fig. ^. 

It seems quite natural to assume that the states be- 
longing to the unitary class in the quasi-lD geometry 
remain to be in the same class as M — > oc. This would 
mean that all the states away from the band center are 
localized. A state at E = 0, if it exists, should not be 
localized in 2D. It follows both from the recent result |2jJ 
that the state at E = is delocalized for odd M under 
open boundary conditions in the y direction and from the 
numerical data in Fig. ^. The derealization of the zero- 
energy state is inferred by requiring that the 2D limit 
(M — > oo) should be independent of the boundary con- 
ditions |33| and of the parity of M. The derealization 
at the band center is a consequence of the particle-hole 
symmetry as in the random hopping model |23| ] and will 
be ruined by the onsite-disorder. 
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FIG. 4. Variance of g as a function of the length of the dis- 
ordered region at E = —0.1. vargr approaches for L shorter 
than mean free path, although invisible in this scale. Inset: 
Variance of g at E = 0. The thin curves are the analytic 
result for the unitary ensemble pi| , where varg — > 1/15 as 
L -> (with M < L). For both figures M = 32. 

As pointed out in j|,|o]], the random-flux model has 
the same symmetry property as the Gade-Wegner model, 
and it is natural to expect that the two models share 
the same critical behavior. In the Gade-Wegner model 
the localization length diverges towards the band cen- 
ter, where the DOS p{E) is also divergent as p(E) ~ 
exp(-c- v /ln|l/£|)/|£'| (c: constant) f||. The charac- 
teristic energy scale below which the singularity of the 
DOS manifests itself is then E c = exp(— c 2 ), which can 
be extremely small depending on c. This may explain 
why no singularity was found in p(E) before ||,|2^,^,Q . 
To find the presumably weak singularity, I computed the 
DOS with high accuracy using the recursive method [ p5[ . 
In this calculation a small imaginary number was added 
to the energy (E — » E + ij) , instead of attaching perfect 
wires. This amounts to averaging p{E) over the energy 
interval of order 7. Figure ^| shows the DOS of a system 
of L = 128000 and M = 64 with 7 = 10~ 2 . The over- 
all shape of the DOS is similar to the one obtained by 
the retraced-path approximation [ ^6[ . Notice, however, 
the tiny peak centered at E = 0. Its height grows with 
smaller 7 and larger M (inset) , which is a clear signature 
of the divergent DOS. To determine the precise form of 
the singularity requires further investigation. It is impor- 
tant to note here that 7 is kept large enough to smear out 
the microscopic structure in p{E) near E = 0. Because 
of the level repulsion and of the particle-hole symmetry, 
p{E) vanishes at E = for even M ]37] . It is expected 
that, in the limit M — ► 00, the dip in the DOS at E = 
disappears and p(E) diverges at E = + , in analogy with 
the ID random-hopping model with even number of sites 
p3[ . The moderate smearing due to 7 helps revealing the 
diverging behavior. 
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FIG. 5. Density of states of a system of M = 64 and 
L = 128000 calculated with 7 = 10~ 2 . Inset: p(E) of a sys- 
tem of M = 128 and L = 64000, 64000, 204800, and 256000 
for 7 = 10~ 2 , 10~ 3 , 10" 4 , and 5 x 10" 5 , respectively. 

The discovery of the divergent DOS at E — estab- 
lishes the connection between the lattice random flux 
model and the Gade-Wegner model. The critical behav- 
ior of the latter model is closely related to the model 
of Dirac fermions with random gauge field which 
has the same particle-hole symmetry. This supports the 
conclusion based on the symmetry argument that a state 
at the band center is the only delocalized state for any 
P (0 < p < 1) in the absence of the on-site disorder. In 
models without the particle-hole symmetry, all the states 
should be localized, in agreement with [p~0^ p~2|JT^ ] . 
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